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Abstract 

We show that the recently derived renormalization group improved Yennie- 
Frautschi-Suura (YFS) exponentiation of soft (k Q — ► 0) gluons in QCD is fully com- 
patible with the usual Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolu- 
tion of the structure functions of hadrons in the respective hadron-hadron hard 
interactions. We show how to implement the YFS exponentiation without double 
or over counting effects already implied by the DGLAP equation. In this way, we ar- 
rive at a theory which allows for the development of realistic, multiple gluon Monte 
Carlo event generators for hard hadron-hadron scattering processes in which the 
DGLAP evolved structure functions are correctly synthesized with the respective 
YFS exponentiated soft gluon effects in a rigorous way. 

Submitted to Phys. Lett. B 



Work supported in part by the US DoE contract DE-FG05-91ER40627, Polish Gov- 
ernment grants KBN 2P30225206 and 2P03B17210 and Polish-French Collaboration 
within IN2P3. 



UTHEP-95-1101 
November 1995 



An important step in the implementation of the recently derived |J renormalization 
group improved soft Yennie-Frautschi-Suura (YFS)[0] gluon exponentiation in QCD is 
the proper synthesization of these new infrared (IR) effects ( where by infrared we do 
mean the gluon energy k a — > limit as opposed to the soft gluon regime analyzed by 
Dokshitzer et a/.|| in which k± > ^, R ~ hadron radius ~ 1/ (200MeV) ~ Aq CD ) 
with the required Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [Q evolution of 
the respective hadron structure functions in computing any particular hard hadron-hadron 
scattering cross section, such as p + p — > ti + X, for example. Indeed, a naive use of the 
results in Ref. would lead to double or over counting some of the same effects, rendering 
the precision of any such calculated cross section questionable at best. In what follows, 
we will show how one properly synthesizes the DGLAP structure function evolution with 
renormalization group improved soft gluon YFS QCD exponentiation derived in Ref. |J • 
In this way, we develope the necessary remaining theoretical apparatus needed for the 
construction of realistic hadron-hadron hard scattering Monte Carlo event generators in 
which finite P± multiple soft gluon effects are taken into account at the level of the hard 
scattering subprocess amplitude while the proper structure function evolution is realized, 
on an event-by-event basis. Such a Monte Carlo event generator will appear elsewhere f| . 

Specifically, let us recall the basic result from Ref. (l| for the YFS exponentiated 
differential cross section for the process q + — > q" + { q ) + n(G), where G is a soft 
gluon, 

da exp = ea;p[SUM IR (QCD)] £~ £ / EI "=i / $Le^ + ^-^ 

xp n (h,...,k n )^^ (1) 



where 



J-JqcD — I -J— J QCD 







-iy-k 



-9(K max -\k\) . (2) 



Here, we have defined 



SUM IR (QCD) = 2a s B QCD + 2a s K£ QCD (3) 

from Ref. [l| with the QCD YFS virtual and real genuinely infrared functions -Bqcd, 
Bqcd, and Sqcd as defined therein Q. The hard gluon residuals j3 n (ki, . . . , k n ) are also 
as defined in Ref. for example. At issue is how one synthesizes a calculation such as 
that in ([!]) with the fundamental structure functions for the incoming partons which are 
themselves solutions of the DGLAP equation and hence in general may have in themselves 
some of the effects in (0) already included? 

Indeed, one might think that the usual formula 

a(pp — > ti + X) = I 2^ Di(xi)Dj(xj)da exp (xiXjs)dxidxj (4) 
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could be used to leading log accuracy when the existing structure functions D^Xi), Dj(xj) 
for partons i,j in p,p, respectively, in the literature are used. Here, the 
respective light-cone momentum fractions. But, this is not even true because the DGLAP 
equation is well known to contain all leading logs from the initial state radiation of gluons 
and all attendant collinear infrared effects. Thus, in using (f|) naively one would double 
count some leading log and collinear infrared effects. 

Specifically, in order to use the existing experience on the structure functions together 
with the YFS exponentiation improvement of perturbation theory given by ([!]), one has 
to identify those effects in the YFS formula which are already contained in the DGLAP 
equation for the structure functions and remove them. This we now do. 

Since we have exponentiated the quantity SUMm(QCD) in (01), we start by isolating 
from it and from all of its attendant parts in ([!]) that collinear infrared contribution which 
is generated by the DGLAP evolution of the structure functions. Viewing the DGLAP 
evolution as the complete LL series from the initial state radiative effects, we see that we 
must do the following to synthesize (0) with the existing structure functions in @: 

• From the initial state part of Sqcd, we must remove that part which corresponds to 
the YFS limit of the DGLAP kernel, Sqcd\dglap- This can be computed using the 
definitions in the original paper of Altarelli and ParisiQ] and we get the result that, 
if we wish to use (||) in (f|), we should replace Sqcd in (ED with Sq CD where 

Sqcd = Sqcd - Sqcd\dglap (5) 

for 



SqCD] DGLAP (*) 



\z) 4(1 - zv) 2 f m z q v z z A 
At: 2 L {k\ + z 2 v 2 m? q ) \ k\ + z 2 v 2 m 2 ) 
6{-z)A{l + zv) 2 f m\v 2 z 2 x i 



(kj_ + z 2 v 2 m 2 ) v k\ + z 2 v 2 m i 



with the identifications zk = k z , k° = v^/2,k = (k z , k±) and we have written the 
RHS of (6) for incoming q and q respectively for definiteness; for the qq' incoming 
case, for example, one would need to substitute m q i for m q in the coefficient of 9(—z) 
on the RHS of (6). The superscript nls on Sqcd on the LHS of (0) denotes that it 
is non-leading in the s-channel relative to DGLAP evolution. Thus, Sq CD has the 
property that its integral over the gluon phase space has no collinear big logarithm 
terms proportional to or 9L ^- , where L = In s/ m 2 — 1 in the case of incoming qq, 
for example. To obtain the result (6), one follows the derivation of the DGLAP kernel 
Pcq(z) in the paper of Altarelli and Parisi in Ref. 0] but one substitutes the YFS 
vertex —igT a (2p + k)^ instead of the complete vertex —igT a ^ ll and one retains the 
quark mass corrections through order m 2 in the respective kinematics. Here, g is the 
QCD coupling constant and r a generate the quark (vector) representation of the QCD 
gauge group; p is the final 4-momentum of the respective radiating quark. We stress 
that Sq CD defined in this way is positive definite. The result (H), which corresponds 
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to the replacement B QCD (K max ) -> B^ CD (K max ) = J ^9(K max - \k\)S^ l CD , then 
avoids the double counting of real gluon radiation from the initial state that is already 
included in the DGLAP evolution of the structure functions. 

• This brings us to the possible double counting of the virtual corrections represented 
by ?R,Bqcd in ([[]). We address this problem as follows. The DGLAP evolution 
equation generates the entire leading log series associated with initial state real and 
virtual gluon radiation. Thus, we should also remove from $IBqcd all leading log 
effects associated with the initial state. Here, we stress that the definition of the big 
log L has to be the same as that used in the DGLAP evolution. In practice, this 
means that, in the formula for SUMm(QCD) in eq.(9) of Ref. for the initial state 
contribution corresponding to the case A = s in the respective sum over channels in 
this formula, we should use the result 

vr 2 1 

Btot(s) = y " 2 (7) 

instead of the result given by eq.(10) in Ref. [0. We denote this new result by 
SUMj^(QCD). In addition, in the computation of the hard gluon residuals j3 n in 
([l|) one should remember again to remove from the real part of the integral over 
the first term in eq.(4) in Ref. |1J (this integral can be found in Refs. ]7)), which 
is the initial state part of JJ-Bqcd; all terms involving either and where 
L = lns/m 2 , — 1 here in the qq incoming case, for example. The generalization to 
the qq' incoming case is immediate. Henceforth, we denote this new BIBqcd with its 
initial state s-channel big logarithms removed as $IBq CD . It is Sq CD an d $IBqcd ^ na ^ 
should now be used in computing j3 n in (JJ); we denote the resulting new hard gluon 
residual as /3™ ls , where it is understood that the various finite order perturbative 
QCD differential cross sections in this definition of are the reduced, hard mass 
factorized cross sections. 

When these steps are followed, we arrive at a new representation of the cross section in 
(|) in which Sqq£), SUM$(QCD) and 5tB$, D are substituted for S QCD , SUM IR {QCD) 
and BIBqcd respectively and we write this new version of the result in as 

da' xp = exp[SUMfA s (QCD)] £~ £ / YlU / ^e +i ^ + ^-^ W+"®n 

x^(k u ...,k n )*^. (8) 

We may now use this result (|8|) in (f|) with the existing DGLAP evolved structure func- 
tions without double counting any effects whatsoever. Moreover, the result (H) is finite 
and has a smooth limit as the initial state light quark masses m q — > 0. 

A related point is that, from the formulas for Sqcd (eq.(5) in the first paper in 
Ref. 0) and for Sqcd\dglap h is clear that the terms we are exponentiating in the factor 
exp{SU Mm(QC 'D)} in (0) are indeed the usual ones that all published calculations have 
treated to a finite order in a s when they combine soft real radiation cross sections with 
virtual gluon radiation corrected ones to cancel the respective infrared singularities at 
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k — ► to the respective order in a s ; here, we treat this cancelation to all orders in a s . 
Referring to Refs. [|S], ^ |It], [II], 0, we see that the value of a s that should be used in 
SUMm(QCD) is a s (fi) where \i is the respective hard renormalization scale~factorization 
scale in the hard process- otherwise, the real gluon emission k — > singularity will not 
cancel the respective virtual k — > singularity: for illustration, we recall the result in 
Ref. || for the parton level Drell-Yan |13| process q + q —> [i + + G for a muon pair 



mass M with the IR gluon regulator mass dig — > , eq.(5.3.17) in Ref. || (here, <7o is the 
respective Born cross section), 

a DY (real) = 2a ^ M h {\n 2 (m 2 G /M 2 ) + 3 \n(m 2 G /M 2 ) + tt 2 } (9) 

37T 

and its corresponding virtual correction, eq.(5.3.18) in Ref. ||, 



a DY (virtual) = ^Mla {~\n 2 (m 2 G /M 2 ) -3ln(m 2 G /M 2 ) - 7 --^ + vr 2 }, 
the two of which sum to the total result, eq.(5.3.19) in Ref. ||, 



(10) 



a DY {real) + a DY {virtual) = — — — cr <^ — > , (11) 

which is independent of the infrared (ko — > 0) gluon regulator mass - the fco — > 
gluon infrared singularities regulated by m G in (|9]) and ( [TO]) are cancelled in (|TT| ) at t/ie 



Ziard scattering coupling a s (M). We refer the reader to Refs. || [10| , to the second paper 
in Refs. |TTJ, and to the second paper in Ref. [O for further explicit examples of this 



choice of a s in the cancelation of IR singularities at k Q — ► in finite order perturbative 
QCD calculations. More recently, the reader may also see the analysis in Ref. |T4| for 
further illustrations of this point. Finally, we recall the fundamental analysis of Altarelli 
and Parisi in Ref. in which the regularized DGLAP kernel, P qq (z), generates the 
differential change in the total qq splitting function given by 



dT qq (z,t) = ^P-P gq (z)dt (12) 
Zty 



and 

1 + z 2 3 



^H^TJ^ + i^ 1 -^)) (13) 
where the +function has the standard definition 

o (1-2)+ Jo (1-^) 1 ' 

for any appropriate test function f(z), Cp = 4/3 is the respective quark color Casimir 
invariant and t = In Q 2 /m 2 — 1 for example. Q is the respective hard scale. We focus 
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here on the 5(1 — z) part of dT qq (z,t) which realizes the cancellation of the k — > real 
and virtual gluon singularities in dT qq (z,t) and which enters into it with the squared 
coupling constant of a s (Q), the coupling constant of the respective hard scale Q, in 
complete analogy with the value of a s which we use in SUM IR (QCD). The correctness 
of the value of a s multiplying the delta function in dT qq (z,t) is proven by agreement 
between the predictions of the DGLAP equation and the Wilson expansion for the deep 
inelastic structure functions (see Ref. j|, [15, for example) and by the agreement of 



these predictions with the available data ||17|| . Thus, the value of a s which we use in 
SU Mir(QC D) is well-founded in the published literature. 

Up to this point in our discussion, we have focussed on how one uses the existing 
DGLAP evolved structure functions in conjunction with the renormalization group im- 
proved YFS theory to compute the rigorous hadron level cross section. It is also possible to 
go one step further, to use the explicit application of the YFS theory to the one-loop and 
single bremsstrahlung correction illustrated in Ref. |18| to exponentiate the IR singulari- 
ties in the DGLAP equation itself. This calculation, which together with its application 
will be reported elswehere [H|, yields an entirely new approach to DGLAP evolution in 
which the splitting functions themselves are actually YFS exponentiated. 

A final important issue that can be addressed is the relationship between our YFS 
exponentiated result ([!]) and the soft gluon resummation theory of Sterman [20] and of 
Catani and Trentadue |21|| . We now turn to this relationship. We note that this latter 
resummation theory has recently been used in Refs. |22], [24|] at the respective leading 
log level to discuss the effects of soft gluons in the p + p — > t + i + X cross section 
normalization. 



Specifically, referring to the notation of Refs. [20 , ^T , [22l 23] , |24| we can identify the 



resummed soft gluon cross section of the theory of Refs. |2(], |21| with the cross section in 
(|) via 



exp[SUM IR (QCD)]^ 



n=0 



W) 4 



d 3 qi d 3 q 2 



e E(1 - z > as) da B /d<$>, 



(15) 



where we note the limit 1 — z — > is the limit studied by Sterman, Catani and Trentadue 
in which the hard scattering invariant mass squared fraction z, defined by s'/s, approaches 
its maximum value 1, and where E(l — z, a s ) is the resummation exponent estimated in 
Ref. [^TJ at the leading-log level and in Ref. []2T] at the next-leading log level and used 
in Refs. p2, |23|, ^] at the leading log level for the normalization of a(tt) at FNAL, <i$ 
is the respective Lorentz invariant phase space considered in Refs. [pOj [2~l[| and d&B is 
the respective Born level cross section. For example, the well-known double integral over 
a s (fc^) in E is easily recovered from the lefthand side of (|T5|) by approximating 



D, 



QCD 



(y) 



D QC d(0) = J i±S QOD (k){l-0(K n 



k)) 
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dz f dkl Cts(kl) CF{l _ m _ z)m ^ _ _ + _ _ _ (lg) 



1 — Z J k]_ 71 

so that for the leading log initial state part of SUM IR (QCD) + Dq C d(0) we recover 
exactly the double integral over a s (/r[) = a s ((l — z)k\) — a s ((l — z)Q 2 ) proposed for E 
in Refs. |20, ^TJ (the • • • in ( |16|) stand for the remainder of -Dqcd(O) |l|). Here, in our 



infinite momentum frame (light-cone) kinematics , we identified (1 — z) max = K max /p\ 
and Q 2 = s. We note that, as the authors in Refs. |20| , j2lfl have pointed-out, in allowing 



a s to run in ( |16|) one resums a certain class of leading logs which in our exact series in 
a s (Q 2 ) in ([!]) are generated by the infinite sum on n in ([!]). We note further that the 
remainder of E corresponds to the next leading log terms generated by the infinite sum 
on n in (|l|) over the (3 n when one makes the approximation 

f^n\ J A J- kj ^ q? g 2 u 



oo 



nll\ n 

ni 



n=0 

(17) 

where we define E nl1 to be the next leading log part of E and to be consistent we have 
defined (3' n to be that part of the (3 n in our exact result (^]) whose leading log content is not 
already resummed by integral over the running a s in (|16|). We see that in our approach, 
the results of Refs. PO, pTJ for E nl1 correspond to the approximation 



{E mi T = ffl ^5( Pl + P2 - qi -q 2 -J2 k M(ki, ■ ■ • , k n )^^-f/da B . 

J j K 3 9l 12 ^g) 

We conclude that, at the parton level, the formulas in Refs. PO, |2T| are entirely contained 
in our formula (0) and that they represent an approximation to our result (0) in the limit 
that 1 — z — * 0. 

One of the important consequences of the identification of the Sterman-Catani-Trentadue 
exponent E in (|15D is that we can address the issue of the presence or absence of 
renormalon |25| behavior associated with the leading log and next-leading log trunca- 
tion of our exact result (0) as it is reperesented by E. Specifically, as the authors in 
Refs. [^U], ^T], ^3J] have noted, as 1 — z — > in (P^, the argument of the running 

a s would appear to approach the Landau pole, creating an ambiguity in the non-leading 
terms in the exponent, for example, of the famous renormalon type ^^if one expands 
a s in terms of this one-loop result- we stress that the two and three loop formulas are 
known and they do not seem to support this pole in a s |^B| • On the other hand, we have 
argued in Ref. 0, by the uncertainty principle, that the truly long wavelength gluons , 
with wavelengths much larger than 1/Aqcd, should decouple from the hard process under 
discussion here. Does this decoupling really happen or not? We look at the exact initial 
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state formula for the sum 



SUM IR (QCD)\ 

initial state 

= {2a s ReB QCD + 2a s B QCD {K max )}\ 



initial state 



C F f d 3 k J2 pi p 2 . 2( M 
a s{k ± ){— 1 r) d(K mSLX - k) 



An 2 J k° p\k p 2 k 

iC F f d 4 k. 1 r 2pi + k 



iCr rd^k 1 , 

4tt 2 J Ti K k 2 -m 2 G + ie } V 

\ 2 

k 2 + 2p l k + ie) 



G -r cc v k 2 + 2pik + ie 
2pi + k \ 2 



(19) 



where we have restored the running ansatz of Refs. |20|, ^j] in their estimate of E. We 
stress that for consistency the same definition of a s has to be used for both soft real 
and soft virtual gluons in order for the infrared singularities to cancel. Carrying out the 
virtual k° integral by standard contour methods in the upper complex k° plane, we get, 
for the residue of the gluon propagator = Res(gluon propagator), the contribution 

9 r r I - 9lL [ — (h 2 \( Pl P2 V 

Z,Oi s ixClDQCD\initial state, Res(gluon propagator) 4tj-2 J faO p p k 

(20) 

so that, when we combine this result with the real emission term in (19), we get finally 
the representation 

SUM IR {QCD)\ imtial state = %{ d ^a s {k\){^- - ^-f 

4vr 2 J k o> Kmax k° Pl k p 2 k 

_C £ ffk v 1 . 

2tt J ti ^ y k 2 -m 2 +ie' 

Res(fermion propagators) " 



2pi + k 2p x + k >! - Tr2 



+ 



k 2 + 2p\k + ie k 2 + 2p\k + ie 



) a s (k'l), 



(21) 



where we have introduced the scale k' 2 for a s in the second term in this last equation 
to take into account that, as this term is not infrared divergent and is dominated by 
gluon momenta 0(yfs/2), k'± should also be of this order for consistency. In the second 
term in (21), only the fermion propagator residues , Res(a),a = fermion propagators, 
enter into the respective sum over residues, as we indicate explicitly. We see in the first 
term in (21) that only gluons with energies exceeding the dummy parameter K max (our 
result (0) is independent of K max and 1 — z max ) are actually involved in SUMjr(QCD) so 
that, if we use the resummation ansatz of Refs. |2"0"|, pT| we do not encounter the regime 



A 2 

(1 — z)s — A.q CD in our result for the LL part of E as we may take (1 — z max ) ^> Q ^ D 
without loss of content or predictive power. This is consistent with the earlier observation 
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that (H]), which is a series in a s (Q 2 ), is consistent with the uncertainty principle and hence 
in fact is insensitive to the behavior of a s at scales O(Aqcd) when Q 2 ^> Aq CD . 

In summary, we have shown how one rigorously synthesizes the powerful results of 
DGLAP evolution and renormalization group improved YFS exponentiation in QCD with- 
out double counting. Explicit Monte Carlo event generator data based on our prescription 
will appear elsewhere f5j. 
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